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ABSTRACT: 

We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of 
the gravitational potential between two spherical bodies of any mass ratio. To accurately 
reproduce the tidal evolution of a system at separations less than five times the radius 
of the larger primary component, the tidal potential due to the presence of a smaller 
secondary component is expanded in terms of Legendre polynomials to arbitrary order 
rather than truncated at leading order as is typically done in studies of well-separated 
system like the Earth and Moon. The equations of tidal evolution including tidal torques, 
the changes in spin rates of the components, and the change in semimajor axis (orbital 
separation) are then derived for binary asteroid systems with circular and equatorial mutual 
orbits. Accounting for higher-order terms in the tidal potential serves to speed up the 
tidal evolution of the system leading to underestimates in the time rates of change of the 
spin rates, semimajor axis, and mean motion in the mutual orbit if such corrections are 
ignored. Special attention is given to the effect of close orbits on the calculation of material 
properties of the components, in terms of the rigidity and tidal dissipation function, based 
on the tidal evolution of the system. It is found that accurate determinations of the physical 
parameters of the system, e.g., densities, sizes, and current separation, are typically more 
important than accounting for higher-order terms in the potential when calculating material 
properties. In the scope of the long-term tidal evolution of the semimajor axis and the 
component spin rates, correcting for close orbits is a small effect, but for an instantaneous 
rate of change in spin rate, semimajor axis, or mean motion, the close-orbit correction can 
be on the order of tens of percent. This work has possible implications for the determination 
of the Roche limit and for spin-state alteration during close flybys. 

Keywords: Gravity - Extended Body Dynamics - Tides - Asteroids - Binary aster- 
oids 
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1. Introduction 



The classic equa tions for tid a. 
seminal papers 



MignardI (11979, 



e.g., 



MacDonaldl (119641): 



evo 



198ll )]. revieA vs [e.g., 



( 19991 )]. and textbooks [e.g., 



u tion in two - body s ystems derived or ut i. 



Goldreichl (Il966h: 



BurnsI (119771) 



Goldreich and Sotei 



Murray and Dermott 



Weidenschilling et al. 



19991); 



Danbyl ( 1l992l )] are based 



( 1l989h 



ized in 



966h: 



Peale 



upon the underlying assumption that the two spherical components in the system are 
separated by several times the radius of the larger primary component. While this 
assumption is valid in planet-satellite systems ^ such as Earth- Moon, Jupiter- Galilean 
satellites, and Saturn-Titan, as well as for Pluto-Charon and the majority of binary 
main-belt asteroids (with 100-km-scale primaries), it is not completely accurate for all 



bi nary asteroids, especial 
by 



those in the near-Earth region. Based upon the compilation 



Walsh and Richardson! (120061 ) of measured and estimated binary asteroid component 
size and semimajor axis parameters, nearly 75% of near-Earth and Mars-crossing binaries 
have inter-compo nent separations between 3 and 5 primary radii. An updated compilation 



of parameters by 



Pravec and Harris! (120071 ) including small main-belt binaries, those with 



primaries less than 10 km in diameter, confirms that 75% of binary systems among 
these three populations have close mutual orbits. In addition, double asteroids, those 



systems with equal- size components that were not counted in the above 



as (69230 ) Hermes (|Margot et al. 



2003; 



Pravec et al 



Antiope (jMerline et al.! . !2000! ; !Michalowski et al.! . !2004! ; !Descamps et al.l. 



2003; 



tallies, such 



Margot et al. 



Frostia, (1089) Tama, (1313) Berna, and (4492) Debussy (!Behrend et al. 



2007!). 



20061), (90) 



(854) 



20061), have 



^There are small natural satellites of the outer planets that orbit very close to their 
primaries, but we must keep in mind that these satellites are part of much more complex 
dynamical systems than simple two-component binaries in addition to having negligible 
masses compared to their primaries. 
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separations within 5 primary radii. The favored formation mechanism for nea r-Earth, Mars - 



cross i ng, and small main-belt bina r ies is rotational fission or mas s shedding (IMargot et al. 



2002 



Richardson and Walsh 



spin-up ( iPravec and Harri: 



of absorbed sunlight (iRubincam 



2006 



Descamps and Marchisl . 



20081) most likely due to YORP 



20071). a torque on the asteroid spiii state due to re-emission 



2000 



Vokrouhlicky and Capekl . 



20021 1 . where the typical 



binaries produced have equatorial mutual o rbits with semimaj or axes between 2 and 4.5 



primary radii and eccentricities below 0.15 ( jWalsh et al. 



20081). Though all binary systems 



in these three populations may not have separations of less than 5 primary radii at present, 
if formed via spin-up, these systems likely have tidally evolved outward from a closer orbit. 



Cornplex g eneralized formulae for tidal e volution are pre se nted b y iKaulal ( 1l964j ) and 



MignardI ( 119801 ) as extensions of the work of 



Darwin 



f ll879a 



18801) that account for 



higher-order terms in the expansion of the tidal potential, though, nearly universally, even 
by Darwin, Kaula, and Mignard themselves, only the leading order is applied in practice 
under the assumption of a distant secondary and the negligibility of higher-order terms. 
To date, the most common application of higher-order expansions of the tidal potential 
is in the Mars-Phobos system where tides on Mars raised by Phobos orbiting at 2.76 
Mars radii are causing the gradual infall of Phobos's orbit. As the separation between 
Mars and Phobos decreases, higher-order terms in the potential expansion must gain 
importance. With this in mind, attempts to understan d the observed secu l ar acc eleration 



of Phobos and the past history of its o r bit date back to 



conti nued with 
with Bills et al. 



Smith and Born 



(119761), iLambeckl ( 119791 ). and 



Redmond and FishI (jl964[ ) and have 



Szetd ( 119831 ). among others, 



(120051 ) presenting the most recent treatment of the subject. 



Because many binaries exist in a regime where traditional assumptions break down, and 
because tidal evolution is most important at small separations, we are motivated to examine 
tidal interactions in close orbits. Here, we expand the gravitational potential between two 
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spherical bodies to arbitrary order as well as allow for a secondary of non-negligible mass. 
We then present the resulting equations for the evolution of the component spin rates and 
the semimajor axis due to the tidal bulges raised on both components when restricted to 
systems with mutual orbits that are both circular and equatorial as suggested for small 
binaries formed via spin-up. The effect of accounting for close orbits is examined and 
compared to the effect of uncertainties in physical parameters of the binary system. 



2. Tidal Potential of Arbitrary Order 

The potential V per unit mass at a point on the surface of the primary body of mass 
Mp, radius -Rp, and uniform density pp due to a secondary of mass Mg, radius Rg, and 
uniform density ps orbiting on a prograde circular path in the equator plane of the primary 
with semimajor axis a measured from the center of mass of the primary is 



V 



where G is the gravitational constant and A is the distance between the center of the 
secondary and the point of interest given by 



1-2 



' cos ^p + ' 



a 



a 



1/2 



(2) 



with iIj measured from the line joining the centers of the primary and secondary [e.g., 



Murray and Dermott 



(1l999l )]. In the spherical polar coordinate system (r, 6, (p) shown 
in Fig. [H with the polar angle 6 measured from the rotation axis of the primary and 
the azimuthal angle (p measured from an arbitrary reference direction fixed in space, the 



separation angle ip between the secondary and the point of interest on the primary is 



5 



cosip = COS 6'p 008 6*8 + sin 6'p sin 6*8 COS (0p — 0s) 



(3) 



For widely separated binary systems where the semimajor axis a is much larger than the 
radius of the primary Rp, the potential is expanded in powers of the small term Rp/a such 
that 



The first term is independent of the position of the point of interest and thus produces no 
force on the primary. The second term provides the force that keeps the mass element at 
the point of interest in a circular orbit about the center of mass of the system. The third 
term is the tidal potential 



that is the focus of past studies of tidal evolution where the the bodies are widely separated 
such as in the Earth-Moon system. However, truncation of the expansion of in (jl]) at 
three terms accurately estimates the true potential in ([1]) only for separations exceeding 
5Rp. For smaller separations, as are often found among binary asteroids, higher orders in 
the expansion of V are necessary. 

The full expansion of the potential in (jlj) may be written concisely as the sum over 
Legendre polynomials P^(cosV^), i.e., zonal harmonic or azimuthally independent surface 
harmonic functions, as 




(4) 



U = -G 



M,Rl 1 



(3 cos^ i^-l) 



(5) 



a3 2 
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where the ^ = 2 term of the expansion of V is the dominant tidal term in (|5]). The full tidal 
potential U including all orders becomes 



a ^-^ \ a I 

e—'i \ / 



(7) 



While we will derive the tidal evolution equations in terms of an arbitrary order Table [T] 
lists the order ^ of the expansion necessary for accurate reproduction of the potential V 
at small separations. At 2i?p, the potential must be expanded to at least £ = 6, requiring 
four additional, but manageable, terms in the expansion. This separation is convenient in 
terms of tidal evolution as it is the contact limit of a binary system with two equal-size 
components and a reasonable initial separation for the onset of tidal evolution in a newly 
formed binary system, regardless of component size, esp ecially for systems formed through 



primary spin-up and mass shedding (I Walsh et al. 



2008 



). Proceeding inward of 2i?p rapidly 



requires an unwieldy number of terms in the expansion (e.g., twice as many additional 
terms are needed at l.S-R^ 



^1- 



3. Roche Limit 



The well-known classical fluid Roche limit is located at a = 2.46i?p ( IChandrasekhar 



19691 ) for equal density components, so that if one considers a secondary just outside the 



fluid limit, one must include the Legendre polynomials of orders £ < 4 in the expansion for 
the potential felt by the primary. For solid, cohesionlesj^ bodies (gravitational aggregates 
or so-called rubble piles) modeled as a dry soil, the Roche limit falls approximately between 



cohesionless material has zero shear strength in the absence of confining pressure. The 
interlocking of the constituent particles under pressure, however, can give the material shear 
strength. 
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1.5i?p and 2Rp flHolsapple and Michel bood boosi : ISharmal . |2009|). The cohesionless Roche 
hmit is based upon a binary system that is not tidally evolving, but the secondary remains 
stressed by its self-gravity, rotation (synchronized to the orb ital period), and the d i ffereri ce 



in gr avity from its near, primary-facing side to its far side. 



Holsapple and Michell (12006 



20081 ) illustrate that the mass ratio of the components has a negligible effect on the Roche 
limit, but one would expect that allowing the secondary to have a more rapid spin or 
allowing for higher-order tidal terms due to its proximity to the primary will increase the 
int ernal stresses a nd push the Roche limit farther from the primary, though, as noted 
by ISharmal (120091 ) . these issues have not been studied in detail. 



With a modest amount of coh esion, the secondary may exist within the stated Roche 



20081 ). For the rough properties of a near-Earth binary of 



limit ( iHolsapple and Michell . 

Pp,s = 2 g/cm^ and Rg = 100 m, a cohesion value of < 100 Pa is enough to hold the 



J. 



secondary together at the surface of the primarjo. For comparison, the surface material 
of comet Tempel 1 excavated by t he Deep Impac t miss ion projectile is estimated to 



have a shear strength of < 65 Pa (jA'Hearn et al. 



10^ Pa (IRichardson et al. 



values up to 250 Pa ( ISchellart 



20051 ) and an effective strength of 



2007): fine -grained terrestrial sand is found to have cohesion 



20001 ). Therefore, it is not unreasonable that in the tidal 
field of the primary, the secondary can stably exist at the very least within the fluid Roche 
limit (even if cohesionless), if not also within the cohesionless Roche limit (with a cohesion 
comparable to comet regolith or sand), justifying our later choice to work to order i = 6 
corresponding to a separation of 2Rp. 



^The cohesion needed to prevent disruption scales as the square of both the density and 
size of the secondary. Thus, for a main-belt binary with a i?s = 10 km, the necessary cohesion 
is of order 10^ Pa, similar to monolithic rock. 
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4. External Potential of Arbitrary Order 

The tidal potential Ui of arbitrary order £ >2 felt by the primary, taken from ([7]), may 
be written concisely as 



Qp 0,p Pi (cos ip) 



where = GMp/Rp is the surface gravity of the primary and 



V 



(9) 



The combination (e.pPe {cosip) is the equilibrium tide height, due to the tidal potential of 
order £, that defines the equipotential surface about a primary that is completely rigid 
(infiexible). Because the mass ratio Mg/Mp < 1 and we assume a > 2Rp, the quantity 
Q,p/Rp ^ 1/8 for all binary systems, and typically Q,p/Rp ^ 1- 

For a body with realistic rigidity, the tidal potential physically deforms the surface 



of the primary by a small distance Xe^pRpSe as a function o 



p ^ 1 and Si is a surface harmonic function. iDarwinI ( 1l879al ) and iLovd (119271 ) lay 



posit ion o n the primar y, where 



the groundwork for showing that, in general, the deformation of a homogeneous density, 
incompressible sphere 



Af,p Rp Si — — hi^p 



El 

9p 



h£,pO,pPe (cost/)) 



(10) 



is given in terms of the displacement Love number h^^p (IMunk and MacDonald 



1960) 



2i+l 



H^p 



2(^-1)1 + 



(2£2+4£+3)Atp 



'11^ 
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introducing fip as the rigidity or shear modulus of the primarjo. For bodies less than 200 



km in radius, as all components of binar y asteroid systems are 



stress due to self-gravity gpR ~ Gp^E? 
structures (i.e., the model proposed by 



flWeidenschilling et al 



:he ri gidity /i dominates the 



19891 ) ■ even for rubble-pile 
Goldreich and Saril (120091 )'). such that the Love 



number <^ 1 for small bodies. With /i^.p and (e^p/Rp small, and noting from ( fTOl) that 
^e,p = hi^pQ^p/ Rp, the assumption of a small deformation factor A^^p is justified. 

Of particular interest is the external potential felt by the secondary now that the 
primary has been deformed. It is this external potential that produces the tidal torque 
that transfers angular momentum through the system. Here, we slightly alter our spherical 
coordinate syste m such that 6 now measures the angle from the axis of symmetry of the 



tidal bulge, as in 



Murray and Dermott 



(1l999l ). such that the surface of the nearly spherical 



primary is now given by 



R 



Rpil + ^ Xi,pPe (cos I 



(12) 



The potential felt at a point external to the primary is the sum of the potential of a 
spherical primary with radius Rp and that of the deform ed shell. However, on. 



Murray and Dermott 



y that due 
( 199^ . will 



to the deformed shell, called the non-central potential by 
contribute to the torque. 

In Fig. O the reciprocal of the distance A between the external point (r, 6, (p) and 
a point on the surface of the primary (r', 6', 0') separated by an angle ip, where r' = R 



•^Darwin 



(Il879al ) realized the correspondence between elastic and viscoelastic media and 
provides a generalized form for the deformation of a viscous spheroid, a function equivalent 
to (fTOj) he calls a, that when applied to an elastic spheroid, in terms of rigidity rather than 
viscosity, is equivalent to the expression found here. 
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from f|T2l) . is 



1 1 °° / /? \ ^ 



(13) 



The use of ^' denotes terms based upon the surface deformation rather than the expansion 
of the distance between the points of interest. The non-central potential (per unit mass 
of the object disturbed by the potential) due to the deformed shell with mass element 

oo 

p^R\ ^ A^/,pPr (cos Q') d (cos 6') d(j)' is 



£'=2 



GppRl Yl E ^^'-P [t) I I ^'^ ^'^ 

\ / «/_9 c_n ^ / J J 



(14) 



£'=2 ^=0 

where the double integral goes over the surface of the primary. The integral of the product 
of two surface harmonics like the L egendre polynomia ls over a surface is zero unless 1 = 1' 



such that for a specific order £ > 2 (IMacRobert 



19671), 



3 . . /P„\'+^ 



2£ 

By defining the more familiar potential Love number 



YVpO,p^p(^) PeicosO). (15) 



3 3 1 

which is of a similar order as hi^p, the non-central potential is written in the form 



R Y^'^ 

Ue^nc = - h,pgpCi,p { j Pe{cos6) (17) 
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such that ?7£,nc at the surface of the primary is simply ke^pUe. Because /i 3> gpR for small 
bodies, the Love number ke,p may be approximated by 



gpPpRp 



2(£- 1) 2£2 + 4£ + 3 /ip £-12£2 + 4£ + 3 /ip ' ' 

Taking the external point to be the position of the secondary orbiting at a distance a from 
the primary, the completqj non-central potential per unit secondary mass due to tides 



raised on the primary is 



/ ^ \ -C^+i) 

Unc = -gp ^ ki,pCi,P [jfj ^) 

The non-central potential drops off quickly with increasing separation as the separation to 
the sixth power for £ = 2 and by an additional square of the separation for each successive 
order. The 6 term in the Legendre polynomial accounts for the angular separation between 
the external point of interest and the tidal bulge of the primary. For the specific location 
of the secondary, we define the angle 5 as the geometric lag angle between the axis of 
symmetry of the tidal bulge and the line connecting the centers of the two components. 

^Here, by complete we mean accounting for all orders £. We have, however, limited the 
result to first order in the Love number /c^ p because terms of order A^^p were ignored in f|T3|) . 
These would have produced higher-order terms in the Love number in the final form of the 
potential in (fT9l) . but because we have argued A^^p and k^^^ are both small quantities, terms 
of second and higher order in the Love number are negligible. 



12 



5. Tidal Dissipation Function Q 



In addition to the rigidity jj,, the response of a homogeneous, incompressible sphere to 
a disturbing potential is characterized by the tidal dissipation function Q defined by 



2-kE* 



(20) 



where E* is the maximum energy st ored in the tidal dis tortion and the integral i s the 



energy dissipated over one cycle [see iGoldreichI (119631 ) or 



Efroimsky and Williams! ( l2009l ) 



for detailed discussions]. This definition is akin to the quality factor in a damped, linear 
oscillator and does not depend on the details of how the energy is dissipated. Friction in 
the response of the body to a tide-raising potential plus the rotation of the body itself (at a 
spin rate u compared to the mean motion n in the mutual orbit about the center of mass 
of the system) lead to misalignment by the geometric lag angle 6. 



The geomet ric lag relates to a phase lag by 



no tation follows 
by iGerstenkornI ( 1l955l ) and 



mpq 



-m 6 sign {u — n), where the impq 



K^aulal (119641). and we have im plicitly assumed a single tidal bulge as done 



independent of the tidal frequencies 



3 



MacDonaldl (1l964l ) by using a single positive geometric lag 6 



J. The tidal dissipation function Q, in turn, relates to 



the phase angle as Q^^^^ = |cot e^. 



Impq I 



I ^Impq I + O ( e^^pg ) 



(lEfroimsky and Williams 



20091) 



^The definition of the phase lag ( iKaulal . 119641 : lEfroimsky and Williamsl . l2009l ). when one 
ignores changes in the periapse and node, is e^mpq = — + q) n — muj] Atimpq, where the 
bracketed term is the tidal frequency and Ati^pq is the positive t ime lag in the response of 



the material to the tidal potential. In the potential expansion by 



Kaulal (1l964j ). only terms 



satisfying i — 2p = m and q = survive for mutual orbits that are circular and equatorial 



such that €£, 



mpq 



—m \u — n\ Atlmpq sign {u — n) = —m 6 sign {u — n), assuming a constant 



time lag and a single (positive) value for geometric lag for all viable combinations of impq. 
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provided energy dissipation is weak {Qempq ^ !)• The absolute value of e^mpq is required 
on physical grounds to ensure that Qempq is positive. Since the tidal dissipation function 
is related to the phase lag, a different Qempq technically applies to each tidal frequency. 
Compared to the dominant order i = 2, where only the impq = 2200 term survives in the 
setup of our problem, 



Qimpq = rnS=-Q^^oo = ^Q-' (21) 

in general, where we define Q = (^2200 such that Qimpq for any tidal frequency is proportional 
to a single value of Q. This simple relation between Qempq and the Q of the dominant tidal 
frequency is a direct result of our assumption of a single geometric lag^ndependent of tidal 
frequency. Such a choice may not be the most realistic physical modeO, but does allow for 
simpler mathematical manipulation. Because Q is necessarily positive regardless of the sign 
of the phase lag, we append sign {u — n) to our forthcoming equations, where the spin rate 
oj relates to the tidally distorted component. If w > ra, the bulge leads; if w < ra, the bulge 
lags behind. 



6. Tidal Torques on the Components 

The force on the secondary due to the distorted primary at order I is — MgVf/^^nc; 
and because we have restricted the problem to a circular, equatorial mutual orbit, the 
tidal bulge remains in the orbit plane, and the sole component of the force is tangential to 



Efroimsky and Williamsl (120091 ) 



''In our model, Q varies inversely with the tidal frequency, 
argue in favor of a rheological model where Q scales to a positive fractional power of the tidal 
frequency (at least for terrestrial planets). It is unclear what rheological model is proper for 
gravitational aggregates like binary asteroids. 
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the mutual orbit. Returning to the notation where ^jJ measures the angle from the axis of 
symmetry of the tidal bulgqj, the force at the location of the secondary is proportional to 
— dPe/dipl^^s pointed in the —ip direction. The value of 6 is taken to be positive as 
stated in the previous section such that, for 6 small, the quantity — dPf/dipl^^g is positive 
and the primary bulge attracts the secondary. For a prograde mutual orbit with Up > n, 
the primary bulge pulls the secondary ahead in the orbit; if Up < n, the primary bulge 
retards the motion of the secondary (see Fig. |3]). The resulting torque vector acting upon 
the orbit of the secondary, which is located at position r with respect to the center of mass 
of the primary, is given by Te^p= r x {—Ms VUe^-ac)- Thus, the torque vector, in general, is 
proportional to — dPi/dipl^^g (^^ x f^. As defined, the direction (sign) of -i^ x f depends 
on whether the tidal bulge leads or lags, and we indicate this in the magnitude of the 
torque via the term sign [u — n) such that the torque on the orbit of the secondary due to 
the £*'^-order deformation of the primary is 



p 

, GM^ f a \ -2(^+1) / SP^cos^, 



Rp \RpJ \ di)p 



(22) 



where 5p is the geometric lag angle between the primary's tidal bulge and the line of centers. 
A positive (negative) torque increases (decreases) the energy of the orbit at a rate Fp n. An 
equal and opposite torque alters the rotational energy of the primary at a rate — FpWp such 
that the total energy E of the system is dissipated over time at a rate E = — Fp {up — n) < 
as heat inside the primary. Though energy is dissipated, angular momentum is conserved 
due to the equal and opposite nature of the torques on the orbit and the rotation of the 



^In this notation, the tidal potential in ([7]) deforms the shape of the component according 
to ( fTOj) and produces the external potential (fTOl) all in terms of the single angle ip- 
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primary. Conservation of angular momentum results in the evolution of the mutual orbit 
and is discussed in the following section. 

A similar torque arises from tides raised on the secondary. By the symmetry of motion 
about the center of mass, the torque g is given by swapping the subscripts p and s in (122|) 
such that 



^ GMl f a ^ dP,{cosij, 



Rs \RsJ \ dip: 

( 

k 



s 



GM^ fR,\ ^^'^ fip f a\ [ dPe (cos ^j, 



Rp \RpJ /is \RpJ \ dip. 



sign (ws - n) (23) 
sign {ujs - n) , 



V's = '5s 



where 6s is the geometric lag angle between the tidal bulge of the secondary and the line of 
centers. This torque changes the orbital energy at a rate Fg n, and the equal and opposite 
torque alters the rotational energy of the secondary at a rate — FgCUg, dissipating energy as 
heat in the secondary at a rate E = — Fg (cjg — n). Torques on the primary and secondary 
weaken for higher orders of i and increasing separations, as expected, and do so in the 
same manner as the non-central potential in ( |T9l) . Once the rotation rate of a component 
synchronizes with the mean motion of the mutual orbit, the associated torque goes to zero 
due to the sign {00 — n) ternj^. Note that we have ignored interactions between the tidal 
bulges of the components as these will depend on the square (or higher powers) of the Love 
numbers, which we have argued are negligible (see Footnote 5). 



^If the mutual orbit were not circular, a radial tide owing to the eccentricity would 
continue to act despite the synchronization of the component spin rate to the mean motion. 
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7. Spin Rate and Semimajor Axis Evolution for Close Orbits 

During tidal evolution, angular momentum is transferred between the spins of the 
components and the mutual orbit. For simplicity, assume that the primary and secondary 
have spin axes parallel to the normal of the mutual orbit plane and rotate in a prograde 
sense. Then, the torque on the distorted primary alters its spin with time at a rate 
u}p — — Fp/Jp, where Ip — a^MpR"^ is the moment of inertia of the primary. The pre-factor 
q; is 2/5 for a uniform density sphere, but can vary with the internal structure of the body, 
and is left as a variable here such that the change in spin rate of the primary is 



ap 1 + K \Rp 



-21+1 



n 



dPe (cos t/'p 



sign (a;p - n) , (24) 



V'p=<5p/ 



recalling that —dPi/dip > for small angles and defining the mass ratio k = Mg/Mp — 
(Ps/Pp) (Rs/Rp)^- Also note that n^, which is proportional to (a/Rp)"^, was introduced via 
Kepler's Third Law, n^a^ — G (Mp + Mg). For rapidly spinning primaries with Up > n, the 
torque will slow the rotation. 

To conserve angular momentum in the system, the change in spin angular momentum, 
given by the torque — F^ p, plus the change in orbital angular momentum must be zero. The 
orbital angular momentum for a circular mutual orbit MpM^/ {Mp + Mg) na^ changes with 
time as (1/2) MpMg/ (Mp + Mg) naa such that conservation requires 



a 

Rr 



2kfr,K , „ 



'2e 



n 



dPi (cos ipp) 



sign {ujp - n) 



(25) 



V'p=5py 



for each order I. For rapidly spinning primaries, the orbit will expand as angular momentum 
is transferred from the spin of the primary to the mutual orbit and, so long as the geometric 
lag remains small, higher orders will cause both more rapid despinning of the primary and 
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faster expansion of the mutual orbit than ^ = 2 alone. A large secondary with k ~ 1 clearly 
causes the most rapid tidal evolution. A small secondary with k ^ 1 will not cause the 
primary to despin appreciably due to the K^-dependence of (Ell), but the separation will 
evolve more readily as (125!) scales 

One can derive the change in the semimajor axis in (1251) by other methods including the 
work done on the orbit and Gauss's formulation of Lagrange's planetary equations. Setting 
the time derivative of the total energy of the orbit —GMpMs/2a, which is GMpMgd/2a^, 
equal to the work do ne on the orb it Ti^pU simplifies to ( 125|) . Using Gauss's formulation for 
spherical bodies (see 



Burns 



(119761 ) for a lucid derivation) and a circular mutual orbit, 



n 



(26) 



where is the tangential component of the disturbing force (per unit mass) from the 



previous section, which is (l/a) dUi^nc/dip sign (wp 



n 



with Ui^nc given by (fTTj) . The 



1 + k) term is not typically present in the Gauss formulation, but is appended here to the 



disturbing functior 



the primary (IRubincam 



10 



due to t 



le non-inertial nature of the coordinate syste m centered on 



19731 ) and is necessary, as stated in 



DarwinI (Il880l ). because the 



primary "must be reduced to rest." One can also argue t he term is necessary to acc ount for 



the reaction of one body to the tidal action of the other (IFerraz-Mello et al. 



20081 ) as the 



^'^Algebraically, from the time rate of change of the orbital energy, d = 2a?E/GMpMs^ 
and the change in orbital energy is further related to the velocity of the secondary r and the 
disturbing force F = — MsVf/nc such that E = i ■ Y = naM^T for a circular mutual orbit. 
Replacing E by naM^T and using Kepler's Third Law, n^a^ = G (Mp + = GMp (1 + k), 
in the expression for a gives ( l26l) for a specific order i. If Ms were ignored in Kepler's Third 
Law, the more familiar form of Gauss's formulation would emerge: a = 2T/n. 
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perturbation is internal to the binary system rather than an exte rnal element (e.g., drag 



force, third body). The (1 + k) term is absent in the formulae of iKaulal (jl964l ). which is 
reasonable if the secondary is of negligible mass, but we wish to allow for an arbitrary mass 
ratio. Substitution of the disturbing force into ( 126|1 produces (125|) . By including the (1 + k) 
term in Kaula's equation (38), evaluating Kaula's F and G functions with zero inclination 
and eccentricity, and recalling that Kaula's eempq = —m 6 sign (wp — n) in our notation, we 
find our evolution of the semimajor axis in fl25|) is a special case of Kaula's generalizatior 
as one would expect. 



we 



The tidal evolution of the secondary follows similarly. The torque on the distorted 
secondary alters its spin with time at a rate Ug = —'^s/Is, 



j sign (ws - n) (27) 

l/'s=<5s / 

^^The product of Kaula's Fimp and Ggpq functions is non-zero for circular, equatorial orbits 
only if i — 2p = m and g = 0. The prefactors on each in the Legendre polynomials listed 
in Table [T] are the values of m for each order i that satisfy i — 2p = m. Thus, the cosine 
terms in the Legendre polynomials we list correspond to impq of 2200, 3110, 3300, 4210, 
4400, 5120, 5310, 5500, 6220, 6410, and 6600. This correspondence allows us to link our 
equations written in terms of Legendre polynomials and a geometric lag to Kaula's equations 
written in terms of the phase lag eempq- While the combinations 2010, 4020, and 6030 satisfy 
i — 2p = m, terms with m = cannot contribute to the tidal evolution of the system because, 
by definition, these terms do not produce a phase lag. These three terms are responsible 
for the ^/'-independent components of the Legendre polynomials with £ = 2,4,6 that vanish 
upon differentiation with respect to ip- 



ki,s 1 /^"^'^ V ^ ^^"^^ 2 [ 9Pi (cOS^s) 
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and alters the semimajor axis of the mutual orbit at a rate of 



(cosies) \ ■ f \ /oc^ 
— ^^^J.gnK-„). (28) 

The Legendre polynomials in Table [1] are written as sums of terms of the form cos mi/) 
where m is an integer (see Footnote 11). Thus, the derivative dPe/dip\^^g is a sum of 
terms of the form sin m6. For small geometric lag angles {Q ^ 1), —dPg/dtpl^^^ > 
and sin m6 ~ m6 such that —dPi/dip\^^g oc Q^^. Because the derivative of a Legendre 
polynomial is proportional to <5~^, only the size ratio of the components and their material 
properties, in terms of their respective fiQ values, determine the relative strength of the 
torques and the relative contributions to the orbit expansion. 








(^£,3 













2e-3 



fJ'sQs 



(29) 



with the relative contribution of the secondary decreasing at higher orders of i and for 
smaller secondaries. Note that the relative strength of the torques is independent of 
the mass and densitjo. For classical i = 2 tides on components with similar material 
properties, the torque due to the distorted secondary is a factor of the size ratio weaker 
than the torque due to the distorted primary. For each higher order in the expansion, the 
relative strength of the torque due to the distorted secondary weakens by the square of the 
size ratio. The changes in the spin rates compare as 



^^However, the absolute strengths of the torques in fl22l) and fl23|) do depend on the masses 
and densities of the components. 
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2e-5 



(30) 



This differs from a generalization of Darwin's result [c.f. iDarwiru ( 1l879bl ). p. 521] because 
we have included the ratio of the Love numbers of the components. At the dominant 
orders, i = 2 and 3, with similar densities, shapes, and material properties, the spin rate 
of the secondary changes faster than the primary. However, interestingly, for i = 4, the 
contributions to the changes in spin rates are equal, and for orders i > 4, the contribution 
to the change in spin rate of the primary is greater than that of the secondary. As with 
the torques, the relative strength of the changes in spin rates weakens by the square of the 
size ratio for each successive order i. For smaller secondaries, the changes in spin rates are 
smaller than for similar mass components, and, for all cases, the process of changing the 
spin of the primary is slower than for the secondary. 

Evaluating the Love number k^^p in ( fT8l) and dPe/dipp from Table [1] explicitly for orders 
i < Q, assuming a small geometric lag angle 6p, and applying f ET]) . the spin of the primary 
changes as 



8 1 n'C'plR'^ 
19 ftp f^pQp 



19/^ 

22 Vi?n 
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Rr 



^ 475 
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Rr 



133 
165 



a 

Rr, 



, (31) 



where n has been replaced with Kepler's Third Law to show the full dependence upon the 
separation of the components a/Rp. Using either ( 1271) or ( l30l) . the spin of the secondary 
changes as 



1 ^'G'PIRI fRs---' 



19 fts /isQs 



""[RpJ \R 



sign [Us-n) 



21 



475 /fty/^y» m/MYAVl (32) 
584 \R^J \rJ 165 \rJ \rJ \' ^ ' 

Assuming similar densities for the components, the change in the spin rate of the primary 
scales as the size ratio of the components to the sixth power (oc k^); the spin rate of the 
secondary scales only as the square of the size ratio at leading order, reinforcing from ( l30|l 
how the spin of the secondary evolves more rapidly than that of the primary, especially for 
small size ratios. 

For close orbits, the separation of the components changes as angular momentum is 
transferred to or from the spins of the components such that the overall change in the 
orbital separation for ^ < 6 is the sum of fl25l) and fl28|) . 



a 
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X 



k{1 + K, 



,1/2 



Rr 



-11/2 



' sign (cjs — n) 




sign [LO^-n] 



sign [10^ -n) 



sign (ws - n) 



sign [i^^-n) 



(33) 



Inside the square brackets, having a secondary of negligible size [R^jR^ — )■ 0) has the same 
effect as having a synchronous secondary [bj^ = n); both make the contribution from the 
secondary vanish. Of course, if one considers the factor outside the square brackets, having 
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a secondary of negligible size makes the mass ratio k negligible, while having a synchronous 
secondary does not directly affect k. The change in the mean motion of the mutual orbit 
follows from Kepler's Third Law and ( l33l) as 



Note that in the above equations (I3T1434I) . any difference in density between the components 
is accounted for in the mass ratio k; otherwise, only the size ratio of the components is 
involved in the terms due to tides raised on the secondary. Obviously, the contribution of 
the secondary is most important when the components are of similar size. Not only is the 
contribution of the secondary weakened because of its smaller size, it should also be despun 
faster than the primary such that its contribution turns off when = n long before the 
primary does the same. Furthermore, each equation has a strong inverse dependence on the 
separation of the components even at £ = 2, and while the inclusion of higher-order terms 
will be strongest at small separations, the orbit of a typical outwardly evolving system will 
expand to a wider separation rapidly. 

8. Effect of Close Orbit Expansion on Tidal Evolution 

Inclusion of higher-order terms for the changes in spin rates and semimajor axis in 
(13TI - I331) speeds up the evolution of the system and decreases the tidal timescales. Using 
up to order £ = 6 compared to £ = 2 results in the spin rates of the components changing 
up to 28% faster at 2_Rp, but falling off quickly with increasing separation (Fig. H]) to less 
than 4% at 5-Rp. The size ratio of the components only affects Ws, where the higher-order 
terms are weaker for smaller secondaries. Similarly, for the change in semimajor axis with 
time, assuming both components are causing the separation to change in the same sense 
(sign (wp — n) and sign {us — n) have the same value), using up to order ^ = 6 (Fig. [5]) 




(34) 
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results in a faster evolution by 21% to 28% at 2Rp and decreases quickly with increasing 
separation. Unlike the changes in spin rates, the largest effect on the evolution of the 
semimajor axis occurs when the size ratio is either unity (equal size) or negligible or when 
the spin of the secondary has synchronized to the mean motion such that the tidal torque 
on the secondary vanishes. The change in semimajor axis with time is least affected by 
the higher-order terms for a size ratio of 0.53 with all other size ratios falling within these 
bounds. According to flMl) for the change in the mean motion with time, the value of h/n 
using higher-order terms compared to £ = 2 has the same form as the change in semimajor 
axis in Fig. [5l 

The strengths of the contributions of the extra terms in the close-orbit correction to 
the change in semimajor axis are listed in Table [21 At 2Rp, higher-order terms with i > 3 
account for nearly 25% of the change in semimajor axis with time. Although the i = 6 term 
is necessary for accurate reproduction of the potential between the bodies to within 1% at 
2Rp, it does not alter the change in semimajor axis with time at the 1% level because of 
the stronger dependence of ( l25l) on separation compared to (E]). The net contribution of the 
higher-order terms in Table [2] decreases by roughly 5% at each value of the separation from 
Tabled] with only the i = 3 term having much consequence beyond 3Rp. 

The total change in the component spin rates as a function of separation, shown in 
Fig. ini is given by integration of the ratio of f l3T|) and f l33|) for the primary and the ratio 
of fl32|) and fl33|) for the secondary. Depending on the size ratio of the components, the total 
change in the spin rate of the primary is enhanced by up to 6% at 2Rp over using £ = 2 
tides only, but not by more than a few percent at larger separations. For the secondary, 
perhaps counter-intuitively, despite the spin of the secondary evolving more rapidly with 
time by adding higher-order terms (Fig. Hj), its evolution with respect to the separation is 
less than when using i = 2 only; the deficit is as large as 22% at 2Rp when the size of the 
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secondary is negligible. This is because for smaller secondaries, the effect of higher-order 
terms on Ug in ( 132|) is reduced, while the effect of higher-order terms on Up is independent 
of the size ratio. Thus, for a rapidly rotating primary, the higher-order terms transfer more 
angular momentum from the spin of the primary to the orbit, expanding the separation 
faster than by £ = 2 tides alone and faster than the spin rate of the secondary changes such 
that the net effect on Aa;s(a) is smaller. 

Integration of (l33|l provides the separation as a function of time. For tidal evolution 
from an initial separation of 2Rp to a final separation of 5-Rp (Fig. [7]), the close-orbit 
correction is strongest at the onset, expanding the separation more rapidly than £ = 2 
tides, but only by about 2% over the same time interval. The contributions from the 
higher-order terms lose strength over time as the separation increases resulting in a net 
effect of expanding the separation by ~1% extra by using i = 6 instead of ^ = 2. From 
Figs. |6] and [71 the integrated effects of the close-orbit correction are small, typically of 
order a few percent; the effects are more noticeable in the instantaneous rates of change 
of the spin rates, separation, or mean motion due to the rapid fall-off in strength of the 
higher-order terms with increasing separation and how rapidly the system tidally evolves 
from small separations. 

Rearrangement and integration of (l33l) allows one to calculate the combination of the 
material properties of the components fiQ (assuming /ipQp = fisQs) and the age of the 
binary At based on measurable system parameters. For brevity, we retain only terms due 
to tides raised on the primary giving 
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with X = a/Rp. Because both terms on the left-hand side of fl35l) are unknown, on e 



may either estimate the material pro perties by assuming binary ages (IMargot et al. 



2002 



2003 



Taylor and Margot 



propertie s (IWalsh and Richardsonl. 



avenues (IMarchis et al 



2008a 



b 



20071). estimate binary ages by a ssuming material 



2006 



Goldreich and Sari 



Taylor and Margot 



20091 ) . or consider both 



20101) • Furthermore, precisely because 



both terms are unknown, assuming a value for one has an intimate effect on the calculation 
of the other as changing one's value by an order of magnitude changes the result of the 
other by an order of magnitude. Thus, when one wishes to find fiQ, for instance, choosing 
an age for the binary injects a great source of uncertainty into the calculation. 

The close-orbit correction enhances the rate at which the separation changes such 
that, to provide the same tidal evolution over the same timescale At, the product fiQ must 
increase to compensate for the inclusion of the higher-order terms. For classical i = 2 tides, 
the denominator of the integrand in f p5|) vanishes such that the effect of including terms up 
to ^ = 6 alters fiQ according to 



^^22-^ ^459 ^584-^ ^165 

and is shown as a function of the final separation in Fig. [HI Note that in Fig. El the 
contribution of the secondary is included in the numerical integration of (!33|) although it 
is not explicitly given in ( 136|) above. Evolution from a close initial separation of IR^ to a 
wide separation of 10i?p results in only a ~1% increase in /xQ over the classical value for all 
size ratios. Thus, the basic 1 = 2 tidal mechanism is sufficient for well-separated binaries. 
On the other hand, if the final separation is smaller, as is the case for most near-Earth 
binaries, the correction is larger, increasing to 5% for evolution from 2i?p to 5-Rp and 15% 
for evolution from 2i?p to 3i?p. When making a coarse estimate of the material properties 
of the system, taking the close orbit into account is not of paramount importance; classical 
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tides will easily provide an order-of-magnitude estimate of fiQ for even the closest of 
binary asteroids, though the result will be slightly underestimated. Complement arily, if 
higher-order terms are included and /iQ is held fixed, the age of the binary must decrease 
by the same factor as in Fig. |8] meaning that i = 2 tides provide an upper bound on ages 
for systems with a given fiQ value. 

The use of higher-order terms up to £ = 6 is sufficient for exploring the tidal evolution 
of binary systems with separations greater than 2Rp. Additional terms with £ > 6 make 
inconsequential changes to tidal evolution at these separations as illustrated by the rapid 
fall-off of the contributions of the higher-order terms beyond 2Rp in Table O Moreover, 
terms with i > 6 leave Figs. HHH] unchanged, only having an effect within 2Rp. Thus, if 
one wishes to proceed inward of 2Rp, simply using orders of up to £ = 6 is insufficient as 
higher-order terms gain importance the closer one proceeds to the primary. Though we 
stated earlier that the number of terms required can rapidly become unwieldy, one can 
approximate their strength. For an arbitrary order i > 2, the term within the square 
brackets of (l33l) is approximately 

allowing additional terms to be included without explicit calculation of the Love numbers 
ke^p or manipulation of the Legendre polynomials. Similar terms follow for the changes in 
spin rates. One must keep in mind that the approximation in fj37|) is only valid so long 
as the small angle approximation holds (sin m6 ~ m6 oc Q^^) with m < i, which requires 
Q > 10 to retain 1% accuracy at m = 6 and larger Q as m increasefl"(e.g.. Q > 20 for 

^■^We have applied ( l2Tl) to estimate the value of Q required. In general, the small angle 
approximation holds to within 1% for Qempq ~ 4 or greater. 
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m = 10). Also, having separations of less than 2i?p requires smaller secondaries, since 
contact occurs at a separation of (1 + Rs/Rp) Rp, which reduces the contribution of the 
secondary due to dependencies upon the size ratio, in addition to demanding consideration 
of the Roche limit for the system (see Section [3]). 



9. Comparison to Measurement Errors 

Take , for examp l e, 66391 (1999 KW4) , the b est-studied of the near-Earth binary 



systems (lOstro et al. 



2006 



Scheeres et al. 



20061 ). Even with exhaustive analysis of 



radar imagery, production of three-dimensional shape models of both components, and 
investigation of the system dynamics, physical parameters of the system are not known 
with extreme precision. The densities of the primary and secondary components are known 
to approximately 12% and 25%, respectively. The uncertainty in the density of the primary 
alone can cause error of more than 30% in Up, Us, and a/Rp according to ( !3TI - [33l) . more 
than the close-orbit correction causes in Figs. H] and [5l The higher estimated density of 
the secondary in the 1999 KW4 system of 2.81 g/cm^, compared to 1.97 g/cm^ for the 
primary, directly affects the mass ratio k applied in the equations of tidal evolution as one 
typically assumes similar densities for the components. Ignoring the density uncertainties, 
this difference in component densities alone causes a 43% change in n that, in turn, affects 
Up by a factor of two and and a/Rp by approximately 40% as well, again, a larger effect 
than the close-orbit correction to tidal evolution. The calculated value of fiQ in ( l35l) is 
affected by density and mass ratio uncertainties in the same way as a/Rp. Furthermore, 
uncertainties in densities and the dependence of the mass ratio k, on density differences 
between the components apply at all separations unlike the close-orbit correction, which 
falls off quickly with increasing separation. 

One must also consider the effect of the initial separation of the components at the 
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onset of tidal evolution, a property that is not known for individ ual systems, but can be 



estimated from simu 



ations of binary formation mechanisms [e.g.. 



Walsh and Richardson 



(120061 ) : IWalsh et al.l (120081 )] and given a lower bound by the contact limit at (1 + Rs/Rp) Rp. 
Assuming evolution over the same timescale, if the system had an initial separation Oi 
instead of 2Rp, the effect on fiQ calculated with classical i = 2 tides raised only on the 
primary is 



13/2 



(3^ 



f^Q2 1 - (ai/af)''/' ■ 

For a final separation Of from 3Rp to lOi^p, unless the actual initial separation is within 
10% of the final separation (>0.9af), the value of fiQ is affected by less than a factor of two 
by assuming an initial separation of 2Rp. Using up to £ = 6 and allowing for tides raised on 
the secondary with any size ratio do not cause a significant difference in this result. 

A similar result is found for the dependence on the final separation of the components, 
which one typically takes to be the current separation. If Of is the final (current) separation, 
then changing the separation to Of/ due to, say, a measurement error causes the calculated 
/iQ value for tidal evolution from 2Rp to change as 



jlQfi _ 1 (af/i?p) 



\ 13/2 

(39) 



We find fiQ is affected by less than a factor of two if the final (current) separation is 
known within 10%. From the dependence on the initial and final separations, it is clear 
that the tidal evolution near the final separation dominates over the early evolution where 
the close-orbit correction is necessary. In fact, if instead of calculating fiQ, one considers 
the time taken to tidally evolve to a final separation Of > 4i?p (by assuming a value of 
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fiQ instead of an age), the evolution of the separation from 0.9af to af takes roughly the 
same amount of time as the evolution from Oi < 2Rp to 0.9af. Thus, precisely when the 
close-orbit correction is most prominent is also when the system requires the least amount 
of time to evolve, which causes the mild effect of the close-orbit correction found in Figs. [7] 
and El 

Returning to a concrete example, for the 1999 KW4 system, using the equivalent 
spherical radius of the prim ary shape raodel, the separation of the components a/Rp is 



known to 3% as 3.87±0.12 ( jOstro et al.l . 120061 ) . By ( l39l) . this small uncertainty can result 



in a roughly 20% error in the calculated fiQ, more than twice the effect of the close-orbit 
correction in Fig. [8] at 3.87-Rp. Together with the dependence of the fiQ calculation on the 
density values for the components, the accuracy of measurements of physical parameters 
in the 1999 KW4 system is more important than accounting for the proximity of the 
components to one another. 



10. Discussion 

We have derived the equations of tidal evolution to arbitrary order in the Legendre 
polynomial expansion of the separation between two spherical bodies in a circular and 
equatorial mutual orbit allowing for accurate representation of evolution within five primary 
radii. Equations written in terms of the Love number ki are applicable to any binary 
system, while equations where the Love number has been evaluated have assumed the 
bodies involved have rigidities that dominate their self-gravitational stress (characteristic 
of bodies less than roughly 200 km in radius). Because higher-order terms cause tidal 
evolution to proceed faster, choosing to ignore them produces upper limits on tidal evolution 
timescales and lower limits on material properties in terms of the product of rigidity and 
the tidal dissipation function. However, we have shown that the correction for close orbits 
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has only a minor integrated effect on outward tidal evolution and the calculation of material 
properties, comparable to or less than the effect of uncertainties in measurable properties 
such as density, mass ratio, and semimajor axis (scaled to the radius of the primary 
component). In the case of outward evolution, the binary system evolves rapidly through 
the range of separations where the close-orbit correction is strongest, so one can safely 
ignore the correction to obtain order-of-magnitude estimates of timescales and material 
properties using the classical equations for tidal evolution. Accounting for higher orders is 
more applicable to studying, famously in the case of Phobos, observed secular accelerations 
and the infall of a secondary to the surface of its primary where the higher-order terms 
instead gain strength. 

Though we have presented the expansion of the gravitational potential and the 
resulting equations of tidal evolution in the context of two asteroids in mutual orbit, the 
essence of this work could be generalized for use in the determination of the Roche limit 
and the study of close flybys. The use of a higher-order expansion of the gravitational 
potential in terms of Legendre polynomials is warranted whenever the separation of two 
bodies is within five times the radius of one of the bodieo (see Tabled]). Historically, in 



body 



2008 



fSridhar and Tremai 


ne. 


1992 


Sharma et al.. 


2006; 


Sharma 



Richardson et al. 



1998 



Holsapple and Michel 



2006 



20091 ). stresses are only considered in the much smaller 



secondary while the primary is assumed to be rigid. For small secondaries, the cohesionless 
Roche limit of 1.5-2i?p is much larger than 5Rs such that higher-order terms in the potential 
expansion are not necessary. However, as larger secondaries are considered {Rs/Rp > 0.1), 



^^The potential felt by the primary requires higher-order terms with £ > 2 if the separation 
is less than 5Rp; the potential felt by the secondary requires higher-order terms with i > 2 



if the separation is less than 5Rs 
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higher-order terms in the gravitational potential will further stress the secondary near the 
Roche limit. Also, with components of increasingly similar size, the assumption of a rigid 
primary is not appropriate; the tidal stress on the primary will deform it from a spherical 
shape and produce an external potential as in Section H] that will in turn further stress 
the secondary. If the components are not spin-locked, tidal torques will also play a role in 
stressing the secondary. Thus, if evaluating the Roche limit for components of similar size 
and/or components that are not spin- locked, one must consider the description presented 



here. For disr uption during a close flyby, or 



passing body (jScheeres 



2001 



Scheeres et al. 



simp 



2000 



y mod ification of the spin state of the 



2004J ). one must consider the proximity 



of the flyby in terms of the expansion of the gravitational potential and whether or not 
tidal bulges can be raised on the components that would produce torques capable of further 
altering the spin state of either component. 

It is also important to remember that the higher-order theory presented here has 
implicitly assumed initially spherical bodies. Extension of this work from spheres to 



ellipsoids or to arbitrary shapes would affect the mutual gravitational potential. 



and angular mome ntum balance, and orbital equations as described by 



ScheeresI t00% 



mear 



and ISharmal (j2010l ) . Once the shape is made nonspherical in the absence of a tidal 
potential, the system is subject to a "direct" torque that naturally occurs from the changing 
gravitational pull felt by the orbiting component due to the nonspherical shape of the other 
component. Accounting for the tidal potential introduces the "indirect" torque described 
here due to the deformation of one component by the gravitational presence of the other 
component. Because the amplitude of the tidal bulge on asteroids, the parameter A in 
this work, can be very small due to its direct dependence on the ratio of self-gravitational 
stress to rigidity, its direct dependence on the mass ratio, and its inverse dependence on the 
separation raised to the third (or higher) power, natural deviations from a spherical shape 



may exceed the amplitude of the tidal bulge. However, one must recall that the direct 
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torques due to a nonspherical shape will change direction as the body rotates under the 
orbiting component tending to cancel the pre- and post-encounter effects of the torque as 
opposed to the indirect torque that is in a consistent direction so long as the bulge always 
leads or lags the orbiting component. It may be important to consider direct torques due 
to natural departures from a spherical shape via the use of shape models: oblate or prolate 
spheroids, triaxial ellipsoids, or vertex models such as those made for the components of 
the 1999 KW4 binary system and other asteroids. 
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Table 1: Order i of Legendre polynomials necessary in the expansion of the gravitational 
potential (E]) of a binary system (with ?/; = 0) to accurately reproduce the full potential 
([T]) to within 1% at separations less than a/R^ ~ 5. If a/R^ is greater than the value 
listed, expansion to the corresponding order ^ suffices. Recall that the fluid Roche limit 
is a/-Rp = 2.46 (see Section [3]). Also note the Legendre polynomials are given in terms of 
cos mip, where m is an integer, rather than the more common form of cos™''?/'. 
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Table 2: Maximum contributions by the successive orders i that aher the semimajor axis 
of the mutual orbit in f l5^ at the separations listed in Table [T] The strengths of the con- 
tributions depend on the size ratio of the components with systems having negligibly small 
secondaries or equal-size components having the strongest contributions from higher order 
terms, which are shown here. Having a synchronous secondary {us = n) also has the same 
effect on d^/d. It is assumed the components have similar /iQ parameters and the effect of 
each component's tides on the semimajor axis are additive. 




Fig. 1. Geometry for the potential felt on the surface of the primary due to the secondary 
orbiting a distance a from the center of mass of the primary. The dashed line is the locus of 
points on the surface of the spherical primary that are separated by the angle ip and distance 
A from the position of the spherical secondary and thus feel the same potential. 




Fig. 2. Geometry for the potential felt at an external point due to the deformation of the 
primary from its initially spherical shape (dashed). Note that here 9 is measured from the 
axis of symmetry of the tidal bulge. 




Fig. 3. When the primary spins faster than the mean motion of the mutual orbit {cUp > n), 
the tidal bulge is carried ahead of the tide-raising secondary. The resulting torques slow the 
rotation of the primary and expand the mutual orbit. When the primary spins slower than 
the mean motion (wp < n), the torques speed up the rotation of the primary and contract 
the mutual orbit. Similar diagrams apply to tides raised on the secondary and whether or 
not ujs > n. Note that ip is measured from the axis of symmetry of the tidal bulge of the 
primary with ip — Sp being the geometric lag angle at the position of the secondary. 
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Fig. 4. Time rates of change of the spin rate of the primary (left) and secondary (right) 
as a function of the separation of the components using all orders up to £ = 6 versus using 
classical i = 2 tides only. The plotted ratio amounts to the bracketed portions of (I3T1) and 
( l32l) for the primary and secondary, respectively. The change in spin rate of the primary 
due to higher-order terms is unaffected by the size ratio of the components. The change 
in spin rate of the secondary is greater for larger size ratios, plotted from top to bottom 
with Rs/Rp = 1, 0.8, 0.6, 0.4, and 0.2. The effect of higher-order terms is always below 30% 
beyond 2Rp and falls below 1% by a separation of 10-Rp for all size ratios. 
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Fig. 5. Time rate of change of the semimajor axis of the mutual orbit as a function of the 
separation of the components using all orders up to £ = 6 versus using classical C = 2 tides. 
The plotted ratio amounts to the bracketed portion of ( l33l) divided by 1 + Rs/Rp with both 
components having similar fiQ parameters and contributing to the evolution in an additive 
sense. The solid curve corresponds to a system with components of equal size {Rg/Rp = 1), 
a secondary of negligible size {Rg/Rp = 0), or a synchronized secondary {us = n). The 
lower bound {dashed curve) is for the size ratio Rs/Rp = 0.53. As in Fig. HI the effect of 
higher-order terms is always below 30% beyond 2Rp and falls below 1% by a separation of 
10i?p for all size ratios. 



1 



cd 

II 



^ 0.9 

CO 

II 

< 



0.8 



1 1 1 r 



1 1 1 r 



1 1 1 r 



Primary 



Secondary 




Fig. 6. Total change in spin rate of the components based on tidal evolution from an 
initial separation of 2Rp using all orders up to £ = 6 versus using classical i = 2 tides. The 
coordinate on the x-axis is the final separation of the tidal evolution. With both components 
contributing in an additive sense in (|33l) . the spin rate of the primary is affected more rapidly 
than in the classical case, while the secondary is affected less rapidly. The maximum change 
in the spin rate of the primary occurs for Rs/Rp = 0.53 {upper solid curve) and the minimum 
is the dashed line at 1 for Rs/Rp = 0, 1 or a synchronized secondary. The lower solid curve 
corresponds to the change in spin rate of the secondary for R^/Rp — )■ 0. For larger size ratios, 
the curve for the secondary moves toward the dashed line at 1. 
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Fig. 7. Evolution of the semimajor axis with time using all orders up to £ = 6 versus 
using classical i = 2 tides. Time is plotted logarithmically and scaled to the time necessary 
for a system to evolve from 2Rp to 5i?p via i = 2 tides. As in Fig. |5l the solid curve 
corresponds to Rg/Rp = 0, 1 or a synchronized secondary and the dashed curve corresponds 
to Rs/Rp = 0.53. Using up to £ = 6 gives a correction of order 1% to classical tides at any 
point in the evolution from 2Rp to 5-Rp. 
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Fig. 8. In terms of the ratio fiQe=e/ fJ'Qe=2 for tidal evolution from 2Rp with both com- 
ponents contributing, systems with Rg/Rp = 0,1 or a synchronized secondary share the 
upper curve and are most affected by the close-orbit correction; a system with Rs/Rp = 0.53 
{dashed curve) is the least affected by the close-orbit correction. Overall, the close-orbit 
correction is roughly 25% at 2Rp and quickly falls off to 5% for evolution from 2Rp to 5Rp 
and to 1% for evolution to 10-Rp for all size ratios. The components are assumed to have 
similar nQ parameters and contribute in an additive sense in (133|) . 



